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Abstract. The basics of the thermodynamic Bethe ansatz equation are given. The 
simplest case is repulsive delta function bosons, the thermodynamic equation contains 
only one unknown function. We also treat the XXX model with spin 1/2 and the XXZ 
model and the XYZ model. This method is very useful for the investigation of the low 
temperature thermodynamics of solvable systems. 



1. Introduction 



Thermodynamic Bethe ansatz equations were first introduced by Yang and Yang 
for the model of repulsive delta-function bosons. Later this method was extended 
to the other Bctho ansatz soluble models such as the XXX model, the XXZ 
model, the XYZ model, delta-function fermions, the Hubbard model, the t-J 
model and so on. In these theories the of strings and holes play an essential role. 
The number of kinds of strings is generally infinite except in some special cases. 
We can construct a set of non-linear integral equations for these models. These 
are very useful for the investigation of the low-temperature thermodynamics of 
soluble models. Nowadays thermodynamic Bethe ansatz equations have been 
derived for almost all solvable one-dimensional quantum models. 

In chapter 2, Yang and Yang's theory for bosons is introduced. Here the 
number of unknown functions is 1. In chapter 3 we treat the XXX model. Here 
the number of unknown functions is infinite. But we can treat this set of integral 
equations analytically in some special cases such as the infinite temperature 
limit, the zero temperature limit, and the high magnetic field limit. In chapter 
3 we treat the XXZ model. The case \A\ > 1 is almost identical to the XXX 
model {A = 1). On the other hand the case |zi| < 1 is more complicated, 
the number of kinds of strings is dependent on a parameter. If (cos""'^ Z\)/7r is 
a rational number, the number of kinds of strings is finite and thermodynamic 
Bethe ansatz equations contain a finite number of unknown functions. In chapter 
5 we treat the XYZ model in zero magnetic field. Here, the number of unknown 
functions is also dependent on parameters. The equations are obtained by a 
slight modification of XXZ model at |Zi| < 1. In chapter 6 we review numerical 
calculations of these Bethe ansatz equations. 



2. Repulsive delta- function bosons 



1 Bethe ansatz equations and uniqueness of the solution 

Here we consider the system 

^ = -Ea^ + 2c^5(x,-x,), (1) 

i=l ' i<j 

the problem of bosons interacting via a repulsive delta function potential with 
periodic boundary conditions. Lieb and Liniger showed that this model is solv- 
able by Bethe ansatz and calculated the ground state energy and elementary 
excitations (Lieb and Liniger (1963), Lieb (1963)). The wave function is as- 
sumed to be a linear combination of A^! plane waves. One must determine A''! 
coefficients A{P) with eigenvalue E at xi < X2 < .•• < 

/ = ^ A{P) exp[i(A;pia;i + kp2X2 + ■■■ + kpNXw)], 
p 

A{P) = Ce{P)l[{kpj - kpk + ic), 

E = Y^kl (2) 

Here P's are permutations and e(P) is the parity of permutation P. The total 
momentum is given by 

N 

/cj's arc called the quasi momenta. If wc apply periodic boundary conditions, 
kj 's should be solutions of the following coupled equations 



t — 1 



The logarithm of these equations is 



fc,L = 27r7,-^2tan-i(^i-^). (4) 
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Here the J^'s form a set of distinct but otherwise arbitrary integers (half-odd 
integers) for odd(even) iV. The quasi-momenta kj are determined by eq. (4) if a 
set of quantum numbers Ij is given. The ground state energy is obtained from a 

Frcdholm type integral equation (Liob and Linigcr (1963), Lieb (1963)). For this 
model Yang and Yang found a non-linear integral equation which gives the free 
energy at a given temperature (Yang and Yang (1969)). This result is important 
because it forms the basic theory for thermodynamics of other solvable models. 

At first we should show that the solution of equation (4) is unique for a given 
set of Ij 's. The equations (4) are equivalent to the following equations 

dB{ki,k2,...,kN) _ Q 
dkj 

B{ki,k2, kN) = ^Yl ~ ^^^^^j + H - 

ei{x) = xtan-\x/c) - |ln(l + (x/cf). (5) 

The extremal point of the function B with N variables is determined by the 
solution of (4). It gives an eigenstate of the Hamiltonian. Let us consider the 
following N X N matrix 

__d^B__ ^ 2c ^ 2c 

~ dkjdkl ~ + c2 + [kj - km? ' C2 + {kj -hf 



This matrix is always positive definite because 

2 

C2 + {kj - kl) 



Ij l<3 



for an arbitrary real vector {uj}. The function iJ is a concave function in N 
dimensional space. Thus a solution of equation (4) for a given set of {Ij} is 
unique. 



2 Holes of quasi momenta and their distribution function 

In the limit of c ^ oo this set of eigenstates is complete. Let us consider the 
following function 

/j(yfc) = fc + |^tan-i^^. (8) 
j 

The position of a hole is defined by 

277 

h{kh) = X unoccupied (half — odd) integer. 

We define the distribution function of holes p'^{k) such that the number of holes 
between k and k + dkis p'^{k)Ldk and the density of particles p{k) such that the 
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number of particles between k and k + dk is p{k)Ldk. In the thermodynamic 
hmit we have 

2-K j p{t) + p^{t)dt = h{k) = fc + 2 y tan-i ^^^p{k')dk'. (9) 

Differentiating this yields 

2n{pik) + p\k)) = 1 + 2 ^.fP^. (10) 

Prom equation (2) one obtains the energy and the particle number per unit 
length 

k^p{k)dk, n= I p{k)dk. (11) 



The entropy of the distribution between k and k + dk is the logarithm of number 
of orderings of Lp{k)dk particles and Lp'''{k)dk holes 

-p{k) Inp(fc) - p'^ik) lnp^{k)]. 

Here we use the Stirling formula ln(n!) ~ n(lnn — 1). Then the entropy per unit 
length is 

/oo 
{p{k) + p^{k)) \n{p{k) + p^{k)) - p{k) Inp(fc) - p''(fc) \up^{k)dk. (12) 
-oo 



3 Thermodynamic equilibrium 

The free energy per unit length / = e — Ts must be minimized under the 
condition that n is constant. / and n arc fiinctionals of p{k) and p'^{k). Thus wc 
should determine p{k) and p'^{k) to minimize f — An. Next we look for parameter 
A such that n takes on its required value. In variational calculus A is called a 
Lagrange multiplier. At the minimum point, the variation of / — An must be 
zero for any infinitesimal variation of functions 

/oo 
(fc2 - A)p{k) - T{{p{k) + p\k)) \n{p{k) + p\k)) 
-OO 

-p{k) Inp(fc) - /(fc) lnp^{k)}dk 

(13) 



One should note that 6p{k) and 5p'^{k) arc not independent. From (10) we have 
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Substituting this into (13), one obtains 

/oo 
Sp{k) 
-oo 

This equation must hold for any arbitrary infinitesimal change of the function 
p{k). Thus 

r/(^)^_,.2 . T r c , p{q) 



must be satisfied at thermodynamic equiUbrium. If we put 
e{k) = Tln{p^{k)/p{k)), this becomes 

e(fc) = fc2 _ ^ _ ly"^ dg__£__ + exp(-^))dg. (16) 

e(fc) can be determined by iteration. This function has the physical meaning 
as the excitation energy for an elementary excitation, as will be shown in next 
section. From equation (10) we have 

From this equation wc can determine p{k). 

If equation (16) is differentiated with respect to the chemical potential A, 
one obtains a linear integral equation 

OA ' + ; ''^c2 + (fc-g)2i + exp(f)) dA ■ ^'^^ 

Comparing this with (17) one finds that 

^fe^ = -27rp(fc)(l + exp(^)) = -2n{p{k) + p\k)). (19) 

The entropy density (12) can be written 

s = J[{p + p^) ln(l + exp(-e(fc)/T)) + p{k)e{k)/T]dk. 

The free energy density is 

f = J{k'- <k))p{k) - T{p{k) + p\k)) ln(l + exp(-^))dfc. (20) 
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Substituting (16) into (20) one obtains 

-T{p{k) + p\k)) ln(l + exp(-^))d/e 
= j Ap{k)dk-T J ln(l+exp(-^)) 

X [p{k) + p\k) - i I ^-^—^p{q)dq]dk. (21) 
Substituting (10) we have a very simple expression for the free energy density 

f = An-Tj^ ln(l + eM-'-^-))- (22) 
The thermodynamic potential density g = f — An is 

giT,A) = -Tj ln(l +exp(-^))g. (23) 

For a given temperature T and chemical potential A, e{k) is determined through 

the non-linear integral equation (16) and obtain thermodynamic potential g from 
(23). All thermodynamic quantities are derived from g{T,A) through thermo- 
dynamic relations 

n=-||, s = -^, e = g + An + Ts,.... (24) 
The pressure pis —g. 

4 Elementciry excitations 

We consider the change of energy and momentum when one particle is moved 
from the thermodynamic equilibrium. Assume that the l-th particle is removed 
and one particle is added between the m-th and m -|- 1-th particle 

{7i >l2>h> ... > In} 

{h > h> ... > Im> I'l > Im+i...Ii-i > Ii+i > ... > In} 
The Bethe ansatz equations for the original state are 

N 



Lfc,=27r/,-^2tan-i^^^ 3 = 1,...,N. (25) 

n=l 



(26) 
(27) 
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For the excited state they are 

Lk'j = 2TrIj - 2 tan-i tiJLti - ^ 2 tan"^ ~ ^" , j ^ / 

c c 

Lk'i = 27r/;' - J2 2 tan-' MjL^. 
The change of momentum and energy are 

AK = kl-ki+J2ikn-kn), 

AE = kl'-kf + J2{k':-kl) (28) 
Subtracting (25) from (26) we have 

^^^■^ + E e2 + (^'-fcO^ (^^^--^fe') 

= 2 tan-i - 2 tan"! . (29) 

c c 

Then equation for the back-flow J(fc) = p{k)AkL is 

(1 + ^) j(.; fc;, ..) - / ^(^.^;_,,)./ (^'; fc.)^^' 

= l(tan-i ^— ^ - tan ^— ^). (30) 

TT C C 

Energy and momentum changes are 

AK{k'i, hi) = k'i-ki+ J J{k; k'l, ki)dk, 

AE{k'i, ki) = k[^ -kf + j 2kJ{k; k[, ki)dk. (31) 

It is clear that AK{ki,ki) = AE{ki,ki) = 0. The differentiation of these with 
respect to k'l is 



dAE 



dk\ 



E r 

^=2.; + y2M^;*i)*, (32) 
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where 

d 

u{k;k'i) = -^J{k;k'i,ki), 

(1 + r}{k))u{k; k'l) - j a{k - k')u{k'; k'i)dk' = a{k - k'l). (33) 
Using this integral equation we find that u{k; k'l) is given by the infinite series 

Substituting this into (32) we get an expression of and as an infinite 
series 

dAK 



+ 1 jdk,dk,a{k[ - k,)^^^^a{k, - fe)3^p^ + ... 

dAE 



dk. 



E f 1 

- = 2k'i + / dkMH - h) 2fci 



+ [ [ dkidk2a{k'i - ki) ^ a{ki - fc2) , , ^ . 2fc2 
J J 1 + V{ki) 'l + r]{k2) 



+ ... 



Prom these infinite series Unear integral equations for and are ob- 
tained, 



dAKm^ki) ,^f,,, , , l_ 



1 dAK{ki,ki) 



■ rj{ki) dki 
dk[ ' y ^ ' "l + r]{kx) dki 



(35) 



On the other hand we have equations for 2-jT{p{k) + p^{k)) and de{k)/dk from 
(10) and (16) 

The solution of non-singular linear integral equation is unique. Thus one obtains 

^^^iiM = 2.(p(fcO + /(M)), 
dk'i dki ' 
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Integrating these we have 



AKik'i, ki) =2Tr [ ' p{k) + p''{k)dk, 



ki 

AE{kiki) = e{ki)-e{ki). (38) 

Thus e(fc) — T \n{p^ (k) / p{k)) has the physical meaning as the energy of the 
elementary excitation. 

5 Some special limits 

5.1 c — oo limit 

In the thermodynamic limit the system is equivalent to ideal spinless fermions 

_ J_ f k^dk 
^~ 2-K J exp(fc2 - A)/T + 1 ~ 
T 

27r . 

by the partial differentiation with respect to k. In this limit c/ (c^ + (fc — (j)^) in the 
integrand of (16) is zero. Then we have have very simple solution e(fc) — k'^ — A. 
Prom equation (17) we have 

Prom (23) the thermodynamic potential per unit length is 

g = ~^ [ dHn(l+exp(-^^^)). (41) 



■ / ln(l + exp[{A - k'^)/T])dk = lim -G/L, (39) 

J L^oc 



2tt _ 

This is equivalent to (39). 
5.2 c= 0+ limit 

In this limit integration kernel c/(c^ + {k — q)"^) can be replaced by ■7rS{k — q). 
Then (16) becomes e{k) = k'^ - A - Tln(l + exp(-e(A;)/T)). Then we obtain 

e(fc)=Tln(exp((A;2-A)/T)-l). (42) 

The Gibbs free energy per site g, pressure p, p{k) and p^{k) are 

ni^(h\ — 

27r' 



5 = -p = T|ln(l-exp(-^^))g, p\k) = ^ 



^^^^ = ^ exp((fc^-A)/T)-l - 



This result coincide with that for ideal bosons. 
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5.3 T = 0+ limit 

Generally speaking e(fc) is a monotonically increasing function of k^. At T = 0+ 
we assume that e{±qo) = 0. Thus we have 

p{k) = for fc^ > Qo, p^{k) = for k^ < ql. 

Equations (10) and (16) for k < qo are 



p{q)dq 



(•qa 



2tt p{k) = 1 + 2c / 

The first equation is equivalent with Lieb-Liniger equation(Lieb and Liniger 
(1963)), 

The above theory was introduced by Yang- Yang (Yang and Yang (1969)). CP. 
Yang solved numerically this equation(C.P. Yang 1970). Very surprisingly it 
seems that this simple non-linear equation gives the exact free energy in the 
thermodynamic limit of the ID repulsive Bosons. The next problem is to find 
thermodynamic Bethe ansatz equation for other soluble models. As the next 
simplest case we treat the 5=1/2 XXX chain. 



3. Thermodynamics of the XXX chain 



6 Wave functions of the XXX chain 

The Heisenberg model was the first model to be treated by the method of Bethe 
ansatz (Bcthc (1931)). In the beginning of the 1930's only the ferromagnetic 
case was considered 

H = -jf2srsf^, + sfsi, + sfsui - 2hJ2s!, 

1=1 1=1 
h>0, Sjv+i = Si. (46) 

This Hamiltonian is defined on a 2^ dimensional vector space. The space is 
classified by the total =J2 ■ The ground state is the state where all spins 
are up and = N/2 

W|0 >= £;o|0 >, Eo = -JN/A - Nh. (47) 
Write a general state \^ > in terms of a wave function /, 

\^ >=Y,f{ni,n2,...,nM)S-^S-^...S-jO >, (48) 
where 1 < ni < n2 < ... < um < N and 2M < A''. The eigenvalue equation is 

-■^X^(l - <5n,+l,nj+i)|/(«l,-,nj + l,nj+i, ...,nM) 
3 

+f{ni,...,nj,nj+i - l,...,n„)| 

+ [Eo -E + {J + 2h)M - Jj2^n,+i,n,+, }/(ni ,n2, -nnM) = 0. 

(49) 

Next we assume that the wave function is of the following form 

M! M 

f{ni,n2,...,nM) = ^ A(P) exp(i^ fep^rij). (50) 

P 3 = 1 
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Choosing 

A{P) = e(P) JJ{e''^''Pi+kp, ^ ^ _ 2e'fcp'), (51) 

M 

E - Eo = 2hM + J^{l-coskj). (52) 

insures that / satisfies (49). If we put e*'^^ = [xj +i)/{xj —i), the wave function 
and energy are written as follows 

M! M 

fim,n2,...,nM) = E^(^) n^^^^^' (53) 
A{P) = De{P) Hixpj - xpi - 2i), (54) 

3<l 

M 

E-EQ = 2hM + 2jy^ . (55) 

^ + 1 

The called rapidities. The periodic boundary condition 

f{xi,X2,...,XM) = f{X2,X3,...,XM,Xl + N), 

implies the xj satisfy 

(:^)" = n(^^^^^)- ^• = i.--^- (56) 
\Xj—iJ t,^Xi—xi—2iJ 

In terms of the rapidity variable, the periodic boundary conditions take on a 
very simple form. 



7 Hulthen's solution for the antiferromagnet 

In actual magnetic substances, the ferromagnetic case is rare. Usual one- dimen- 
sional magnetic substances are antiferromagnetic (the J < case of equation 
(46)). The logarithm of equation (56) is 

M 

2N tan" ^ = 27r/j + 2 V tan" ^ ^' , (57) 
1=1 ^ 

where Ij is an integer (half-odd integer) for odd(even) N — M. The total momen- 
tum is given by 

i^ = 7r(l-(-l)^)/2-|^/,. 

j 
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For simplicity we set N to be even. One can show that the lowest energy state 
in the subspace of total Sz = N/2 — M is given by 

Ij = {M + l-2j)/2, i = l,2,...,M. (58) 

In the thermodynamic limit xj's distribute from —B to B. From equation (57) 
we have 

tan"^ x = n j p{t)dt + j tan'^ ^-^p{y)dy. (59) 
Differentiating with respect to x yields 

11 1 



'^(^) = — 2— T - / -? ^2—i;piy)dy- (60) 

nx^ + 1 J_B -K {x-yY +A 



The energy and magnetization per site are 
E _\J\ 



(61) 



(62) 



This integral equation can bo solved in the case of infinite B. We define the 
Fourier transform of p{x) as follows 

/oo 
e-''''' p{x)dx. (63) 
-oo 

Using the formula / 7r~^n/(x^+n^) exj£){—ixijj)dx = exp(— n|a;|), one can rewrite 
(60) as follows 

p(a;)(l + e-2|'^l) =e-l'^l. 
Thus p{oj) = 1/(2 cosh w) and 

''^"^^ = ^ /I ^''^^P^'''^'^ = Jsech(^). (64) 
Substituting this into (62) we have 

1 



e ■ 



|J|(ln2- -) = -0.4431471 J|, = 0. (65) 



It can be shown that the case i? = oo is the true ground state at /i = and 
that the magnetization is zero. Thus an analytical result for a one-dimensional 
antiferromagnet is obtained (Hulthen (1938)). The first neighbor correlation is 
derived from this result 

< 5f5?+i >= ^(1 - 41n2) = -0.14771573. (66) 
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An analytic expression for the second neighbor correlation function is also known 
for this model (Takahashi (1977)), 

< >= ^(1 - 161n2 + 9C(3)) = 0.06067977. (67) 

This was calculated from the ground state energy of the Hubbard model. 
Hulthen's solution is the ground state of XXX antiferromagnet at zero mag- 
netic field. Griffiths calculated the magnetization curve of this model (Griffiths 
(1964)). des Cloizeaux and Pearson calculated the elementary excitation away 
from this ground state (des Cloizeaux and Pearson (1962)). 

8 String solution of infinite system 

The elementary excitation for the M = 1 case is 

E-Eo = J{l- cos K) + 2h. (68) 

This is the spin wave excitation of the ferromagnetic Heisenberg model. Bethe 
found that the bound state of spin waves exists (Bethe, 1931), 

Xj = a + i{n + 1 - 2j),j = 1,2, ...,n. (69) 
The energy and momentum of this excitation arc 

2nJ „ , „ 1 
+ 

Thus the dispersion relation is 



E = Eo + ^^ + 2nh, i^=iln(^). (70) 



2J 

E = Eq + —(1 - cosK) + 2nh. (71) 
n 

These string solutions play an essential role in the thermodynamics of soluble 
models. The author obtained the thermodynamic Bethe ansatz equation for the 
XXX model (Takahashi (1971a)). The wave function is 

M , . 

' xpj + ^ N 



/(ni,n2,...,nM) =Xl^(^)n(: 



A{P) = De{P) Hixpj - xpi - 2i). (72) 

3<l 

We consider the wave function in the equation (72) assuming that N is infinity. 

Particle coordinates rij move from — oo to oo. Assume that > '^k^ > ■■■ > 
^kn- The wave function is written as follows 

M M 

/(ni,n2,...,nM) = (^i^2..^m)"^ ^ ^(F) HCll^^)"'^'""'' 

P j=2 l=j 

A{P) = e(P) Hixpi - xpi - 2i), zj = e^'^' = (^). (73) 

3<l ■' 
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In the infinite system / satisfies the following boundary conditions 
I lim f{ni,n2,...,nM)\ < 00, 

I lim f {m, n2, ...,nM)\ <oo, (74) 

From this condition we find that \ziZ2--zm\ = 1 and that A{P) = if one of 
in^lj^pd is greater than 1. Prom the normalizability condition of the wave 
function we have 

^(7)^0, A{P^I) = 0, \ziZ2...Zn\ = 1. (75) 
These conditions are satisfied only if 

Xj = a + {n + 1 -2j)i, j = l,2,...,n. (76) 

Prom \ziZ2--Zm\ = \ a-nl I ~ ^® must be real. The condition 

a — ni 

is automatically satisfied for general n. Thus strings with arbitrary length are 
possible for the XXX Heisenberg chain. Moreover, the following type of string is 
impossible 

a;i = a + (2 + X2 = a + 0i,X3 = a - pi,X4 = a - {2 + (77) 
where a,/3 are real and (3^1. 



9 String hypothesis for a long XXX ring 

We assume that all rapidities Xj belong to bound states with n ~ 1,2, .... For 
a bound states of n-x's the real parts of all x's are the same and imaginary 
parts are (n — (n — 3)i, — (n — l)i within the accuracy of 0{exp{—SN)). 
This assumption seems to be too strong and there arc some counter examples in 
special cases. This is a very controversial point of thermodynamic Bethe-ansatz 
equations for soluble models, except in the repulsive Boson case, which has no 
the string solution. But equations obtained using the string hypothesis seem to 
give the correct free energy and other thermodynamic quantities. 

Consider the case where M„ bound states of n — a;'s exist. We designate a;'s 

as 

xy, a = 1,2,..., Mr,, 
x^J =xl+i{n + l- 2j) + deviation. (78) 

Prom (56), 

e-(x-)= n ^( VfM^^^^) n 

{m,l3)=it{n,a) j'^j 
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j = l,2,...,n. (79) 

Here e{x) = {x+i)/{x — i). The last product is delicate because the numerator or 
denominator may become very small. If we take the product of these n equations, 
these delicate terms are canceled and we have 

n 

e^«/r^) = ^e'^«'')= H Er,m{x", - x^) , (80) 

where 

{^( \n-m\ )^ ( ln-ml+2)^ ( |n-m|+4)---^ ( n+m-2 )^( ri+m ) 
for n / m, (81) 
e2(|)e2(f)...e2(5^)e(|j) for n = m. 

The logarithm of these equations gives 

N9{xZ/n) = 27,12+ J2 ^nm{xZ-x^), (82) 

e{x) =2tan-\x), (83) 

+ ^^(| n-m|+2 ) + - + 2^( n+^-2 ) + 

for n ^ m, (84) 
2^(1) + 2^(1) + ... + 2e(^)+e(^) for n = m. 

is an integer (half-odd integer) if A'' — is odd (even) and should satisfy 

^ 00 

\I2\<2^N-l-Y,tnmMm), tnm ^ 2Min{n,m) - 6nm- (85) 
m=l 

We can prove the number of sets is — C'm-I' under the condition M = 
Er=i "^^n- Here C£ is the binomial coefficient defined by N\/{M\{N - M)\). 
For details, see Appendix A of (Takahashi (1971a)). The energy of this state is 
given by 

Em}) = N{-h -'-) + E( fan)2 + ^2 + (86) 

We can construct wave functions through (48), (50) and (54) a,t S = Sz = N/2~ 
M. The wave functions for Sz = S—1, S — 2, —S are obtained by applying the 
operator S^^^^^^^. The energy for these cases are E—2h, E—Ah, E~2{N -2M)h. 
Then the total number of states which are generated by the string assumption 
and the descending operator at Sz = N/2 — M is C^. Therefore the total number 
of states is 2^. This coincides with the true total number of states. It is expected 
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that all eigen functions constructed in the above way should be a complete set. 
The partition function of this system is written as follows 

l-exp(-2VT) 2.-P[-^ (87) 

M=0 ' ' {/J} 

The free energy is given by G = — Tln^. We define functions hn[x) by 

K{X) = Onix) - N-^ J2 ^^rnix - C)- (88) 
(m,a) 

We can define holes in n-string sea by the solution of 

27rJ^/7V = (89) 
where Jg are omitted integers or half-odd integers in the region. 

10 Thermodynamic Bethe-ansatz equations for XXX 
chain 

In the thermodynamic limit, define distribution functions of n-strings and holes 
of n-string as pn{x) and pl\{x). The number of strings and holes between x and 
a; + da; is pn{x)Ndx and p'^{x)Ndx, respectively. Thus we have 

/X /"OO 
Pn{t) + pl{t)dt = 6n{x) - Onm{x-y)Pm{y)dy. (90) 

Differentiating by x, we obtain the integral equation 

/oo 
Tnm{x - y)pm{y)dy, (91) 

where Tnm{x) is a function defined by 

{a\n-m\ix) + 2ain-m\+2{.x) + 2a\n-jn\+4ix) + - 
+2an+m-2{x) + an+m{x) for n ^ m, (92) 
2a2{x) + 2a4{x) + ... + 2a2n-2{x) +a2n{x) for n = m. 

an{x) is a function defined by 

1 n 

a„{x) = — 2' (^o{x) = S{x). (93) 

TT x + n 

The energy per site is 

J ^ POO 

e = -(7 + /i) + ^/ gn{x)pn{x)dx, 

gn{x) = 2nJan{x) + 2nh. (94) 
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The total entropy per site is 



' = E Pr.{x) HI + ^) + p>:,{x) ln(l + 



(95) 



e — Ts should be minimized at the thermodynamic equilibrium. Consider the 
functional variation of the free energy with respect to Pn{k) and Pn{k), 



oo „ 

= 6e-TSs = ^ dx 



[g^{x) - Tln(l + P^)]5pn{x) - Tln(l + ^-^)5pt{x). 



(96) 



From equation (88) we have 

5p';,{x) = -Spn{x) -J2 ^"'"(^ ~ y)Spm{y)dy. (97) 
Substituting these into (96) yields 

= ^E /{^-ln^n(^) + E / Tnm{x-y)Hl+ri;;,\y))dy]6pn{x)dx, 

(98) 

where r}n{x) = p'^{x) / pn{x) . Thus we have integral equations for an infinite 
number of unknown r]n{x) 

Inrinix) = + ^ T„„ * ln(l + ^^-^(^r)). (99) 

In the theory of Bethe ansatz equations we encounter very frequently the inte- 
gration of the foUowing type 

/oo 
a{x - y)b{y)dy. 
-OO 

This is a convolution of two functions a{x) and b{x). In the above and hereafter 
we write the convolution as a * b{x), 

/oo 
a{x - y)b{y)dy. (100) 
-oo 

The free energy per site becomes as follows 

f = e-Ts = -{- + h) + J2 J 9nPn - T[pn ln(l + Vn) + Pn H'^ + Vn^)]dx. 
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We eliminate pjj using (91) 



oo 



+p„[ln77„ - ^ - T„„ * ln(l + v~')]dx. (101) 



7 oo . 

/ = -(t + /l) - / «"(^) ^'^(l + ^n\^))dx. 



On 
T 

Using (99), one sees the inside of the bracket on r.h.s. is zero. So we have 

(102) 

1 

From the n = 1 case of (99) we have 

, /I N 2-KJai(x)+2h v~^, \ , /I -In /-.^ox 

ln(l + 77i) = '-y + 2Z(«'-i + * 1=^(1 + 1i )• (103) 

1=1 

Operating J dxs{x) on this equation yields 

J dxs{x) ln(l + ?7i) = -^ J s{x)ai{x)dx + ^ 

oo „ 

+ / a((a;)ln(l + ?7f i(a;))dx. 
Then equation (102) is transformed as follows 

/ = J(ln2-i)-r / s(x)ln(l +?7i(x))dx. (104) 



Solutions rjn of (99) are functions of x, J, T and h. Differentiating (99) with 
respect to J yields 

2TTJan{x) I drjn ^iT^ 1 drjm 

T r]n dJ (1 + r/r„)?7„j aj 

Comparing this with (91) we have 

^" = 2^(1 + ,„),„ a/' ^" + ^" = 2;^^- (1^^) 

By the definition (92), 

Ol * (T„_i^m + r„+i^m) — (oo + 02) * Tn^m = (<5n-l,m + <5n+l,m)ai- (107) 

Using equations (99), (107) yields 

, N , / N 2nJai(x) , , 

(ao + 02) * In 771(2;) = TfT-^ + «i * 1^(1 + ^2(2;)), (108) 
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(ao + 02) *lnr/„(a;) = ai *ln(l + ry„_i(a;))(l + r/„+i(a;)), n = 2,3, .... (109) 

Equations (108) and(109) are not complete to determine all of the r]n{x), as they 
do not contain h. Take the n = 1 case of (99) 

2nJai(x) + 2h , 1, ■r^, n , /, ix 

lnr?i = +a2*H^ + Vi ) + 22{aj-i + aj+i)*ln{l + r]J^). (110) 

Substituting (108), (109) we can eliminate r]j,j < n for a given integer n, 
2h 

— = an* Inrjn+i - fln+i * ln(l + rjn) - an+2 * ln(l + r?^+i) 

00 

- ^ (ai_i+ai+i)*ln(l + 77;-i). (Ill) 

l=n+2 



Thus we have 



2/1 

ln?7„+i = —+ai* Inrjn + 02 * ln(l + 77~+i) 



+ ^ (a,_„_i +a;_„+i) *ln(l + 77/). (112) 

l=n+2 

For large n, ln(l + r]~^) ~ o(n~^) and therefore: 

lim lnr?„+i - ai * ln77„ = (113) 



or 



lim ^ = f . (114) 



Thus the following equations determine r/„, 



where 



2ttJ 

lnr]i{x) = -^s{x) + s* ln(l + r]2{x)), (115) 

lnr?n(a;) = s* ln(l + r?„_i(a;))(l + 7?„+i(a;)), (116) 

lim ^ = ^, (117) 



1 TTT 

s(a.) = -sech(-). (118) 
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11 Some special cases and expansions 

11.1 J IT case 

In the limit J/T — > and h/T > we can expect that ?7„(x) is independent of 
X, because there are no x dependent terms in equations (115), (116 and (117). 
As / dxs{x) = 1/2, equations (115,116,117) become 

r?2 = (l + r,„_i)(l + ,y„+i), (119) 
r]1 = l + T]2, lim ln??„/n = 2h/T. (120) 

n— >oo 

Equation (119) is a difference equation of second order. It is similar to a differ- 
ential equation of the second order and contains two arbitrary parameters. The 
general solution of this equation is 

Vn = ( r - 1- (121) 

Parameters a and z are determined by (120) and we have a = z,z = exp{h/T) 
and 

,sinh[(n+l)Vr],2 T n , n 
smh[/?./ i J 

77„ = (n + 1)2 - 1 for h = 0. (122) 
Substituting this into (104) we obtain the free energy, magnetization and entropy 

/ = -T ln[2 cosh /i/T], m = 2s^ = -df/dh = tanh/i/T, 

5 = -df/dT = ln[2 cosh(/i/T)] - {h/T) tanh(/i/r)- (123) 

At /i = the entropy per site is In 2. This corresponds to the fact that the 
number of states per site is two. 

11.2 High temperature expansion or small J expansion 

For the XXX chain, we can perform the high temperature expansion of the free 
energy density from the definition 

//T = -N-^ lnTrexp(-H/T). (124) 

This is expanded as a power series of 1/T. Assume that Ti = TLq + Tii where 
T-Lq and Tii commute each other. Then the exponential operator of Ti, can be 
expanded as follows 

eM-n/T) = exp(-7^o/T)(l - T'^^ + T'^^ - +...) . (125) 

Thus, 



where 
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^ <Hl>-3<ni >^< ni>+2< Hi >3 

+ 3!7vr3 ^^''^> 

TreM-no/T)X 
< ^ >= TreM-no/T) ' ^''^^ 

In the beginning of 1930's only the ferromagnetic model was considered. If we 
set 

JV 

no = -2hY,S^, m = -jJ2SiSr+i + SfSl,+S^SU„ (128) 
1=1 1=1 

for the Hamiltonian (46), we obtain the J/T expansion of free energy at fixed 
h/T 

f/T =-\n{2 cosh h/T) - ^tanh^ (h/T) 

(3 + 2tanh2(/i/r) - 3tanh''(Vl')) + 0{{J/Tf). (129) 



32T2 



The calculation to higher orders can be done by the use of linked cluster expan- 
sion. Higher order terms are polynomials of tanh(ft,/T). 

Apparently the expression of the free energy in (123) coincides with the first 
term of the above expansion. Writing ln(77„ + 1) as the expansion 

1 °° 7 

ln(r?„(a;) + 1) = ln[-— r] +Y,fi%f, 

^ 1=1 ^ 
sinh^(Mn + l)/T) 
" ~ smh{hn/T) sinh(/i(n + 2)/T) ' ^ ' 

we obtain an expansion of lnr/„(a;) 

lnr?„(x)=ln-^ + [^K/W 

+ + (a„ - al)^J^) + 0{[in (131) 

Substituting these expansions into (115,116,117,118) and taking first order terms 
in J/T, linear integral equations for fn\x) are obtained 

aif['\x) - s * f^'\x) = 2ns{x), (132) 
ar^fi'\x) - s * ifi'},{x) + = 0, (133) 

lim ^IhM^ = 0. (134) 

n— J-cx) 77, 
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The r.h.s. of these equations are inhomogeneous terms of the integral equations. 
The Fourier transform of these equations are 

(el-l + e-l-l)a„/(i)(a;) = fj\{u;) + f^%{uj). (135) 

The general solution of this difference equation is 

/Ph = 

sinh((n + 2)fe/r) _ sinHnh/T) 

^^^smh{{n + l)h/T) sinh((n + l)/i/r) ^ 



sinh((n + 2)Vr) _ sinh(nVr) 
^ '^sinh{{n+l)h/T) smh{{n+l)h/T) ^ 



From the boundary conditions we have 



(136) 



A(lj) = — ^^riT^, B(uj) = 0. (137) 
^ ' cosh/i/T' ^ ^ ^ ' 



Thus 



f(^)(uj) = TT sinh3Vr _ sinhfe/T 3,^1 
•'1 ^ ' cosh{h/T)^smh2h/T smh2h/T ^' 

Substituting this into (104) we obtain the second term of the J/T expansion 
(129). The higher order terms can be calculated by solving the linear integral 
equations for fn\fn\ •••• The equations are similar to (132,133,134) except for 
the inhomogeneous terms, which are given by lower order fn^ . 



11.3 Low temperature limit 

At low temperature, lnry„ diverges as 1/T. So we should define the following 
functions 

enix) =Tljir]n{x). (139) 
The integral equations become 

ei{x) = 2-kJs{x) + s * Tln(l + exp(^^)), 
er.{x) = s * Tln(l + exp(^Ii^))(l + exp(^I^)), 
lim fliM = 2h. (140) 

n— *-oo n 
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The free energy expression becomes 

J oo . 

f =-{j + h)-TY^ / + exp(-e„(a;)/T)da; 

n=l 

= J(ln2-i)-r j s{x)ln{l+exp{ei{x)/T))dx. (141) 

The T = Hmit of these equations is 

ei{x) = 2ttJs{x) + s * ejt(.T), 
en{x) = s * (e+_i(x) + e++i(x)), 

lim = 2h, (142) 

n— *oo 77, 



(a;)(ia;, 



f=-d + h) + f2[ a„ix)e-ix)dx = J(ln2 - i) - / s(a;)e+ 

n=l '' 

'"^^^-\0, for e„(a;) < 0,. 

e„-(x)^(°V , iJn' (143) 

i en(a;), for e„(a;) < 0. ' 

In the ferromagnetic case J > we have 

e„{x) = e+(a;) = 27rJa„(a;) + 2/in, n = 1, 2, (144) 

and therefore / = —{j + h). This is the ground state energy of the ferromagnetic 
case. 

In the antiferromagnetic case J < we have 

en{x) = e+{x)=an-i*e+{x)+2{n-l)h, n = 2,3,.... (145) 
The equation which determines ei is 

ei(a;) = -27r|J|s(a;) + /i+ / R{x - y)ei{y)dy, ei{±B) = 0. (146) 

J\v\>B 

In the limit oi h ^ B becomes infinite. We have ei{x) = — 27r|J|s(a;) and 
/ = J(ln2-i). 
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11.4 Fugacity expansion 

In the case of very large h the free energy ean be expanded as a power series 
oi z = ex.p{—h/T), z is the called the fugacity. From equation (94) and (99) we 
have expansions of r]~^ as follows 

Vn^ = oxp(-^a„(a:)) exp[-r„™ * {tj-^ - ir/-^ + -...)]. (147) 

The expansion of ryf ^ and r?^^ up to z'^ is 

VT^ = exp(-^ai(a;))(l - ^ 02(2; - y) exp(-^ai(j/))dy) + 0{z^), 

r,2-i = z4exp(- — a2(a:)) + 0(zf^). (148) 

As r]~^ becomes small, (102) is more convenient than (104) 

J °° r 1 

f = -_-h-TY, dxan{x){rj-' - -rj-^ + -...). 

n=l 

Substituting (148), we obtain 

f = - h - z^T j ai{x) exp(-^^ai(x))c;a; 

im ( f / \ / 27rJ , , , 1 , , , An J , , , 
-z T| / 02(0;) exp( — —ai^x)) - -ai(a;)exp( — —ax[x))dx 

- J dx J dyai{x)exp -^^{ai{x) + ai{y)) a2{x - y)^ + 0{z^). 

(149) 

Putting X = tanu, y = tanv, 
{f+^ + h)/T 



^/o(i^) + ^^{-ie-2^/o(2K) + e-^/2jo(if/2) 

^ g2K(l -cos a^i cos ) (1 _ cos ^OS W2 ) 



dijJidijJ2^ 



1 — 2 cos wi cosa;2 + cos^ wi 
+0(^^), (150) 

where if = J/T,u!i ~ tt + u + v,uj2 — u~v and Io{x) is modified Bessel function. 
This result is the same as that of Katsura (Katsura (1965)). 

The strings are stable in a chain of infinite length in the case of very few down 
spins. Some counter examples of string assumption are found in some special 
limit. Nevertheless the thermodynamic Bethe ansatz equation seems to give the 
exact free energy in the case where the density of down-spins is comparable to 
that of up-spins. This non- linear integral equation contains an infinite number of 
unknown functions. To solve this equation one needs to do numerical calculations 
by computer. 



4. Thermodynamics of the XXZ model 



12 Symmetry of the Hamiltonian 

In this section we consider the following Hamiltonian 

N 

n{J, A,h) = -JY^ ^r^r+i + SfSl, + AStSt^^ - 2ft ^ St, 

1=1 1=1 
h>0, Sjv+i^Si. (151) 

This Hamiltonian contains an additional parameter A. The case Z\ = is called 
the XY model, which can be mapped to non-interacting fermions making it 
possible to calculate many physical quantities (Lieb Schultz and Mattis (1961), 
Katsura (1962)). The case Z\ = 1 is the XXX model and was treated in the 
previous chapter. The limit of very large A is the Ising model. The generalization 
of Bethe's method to Z\ ^ 1 was done by Orbach and Walker (Orbach 1958 
Walker 1959). Yang and Yang investigated the ground state of this model in 
detail (Yang and Yang (1966)). Bonner and Fisher investigated this model using 
the diagonalization method up to A'' = 12 (Bonner and Fisher (1964)). In this 
Hamiltonian the magnetic field is applied in the z-direction. For a magnetic field 
in a different direction, the exact solution is not known. 
Let us consider the following unitary transformation: 

oo 

n{J, A, h) = UiH{J, A, -h)U^\ ^1 = n = ^i"^- (^^2) 

1=1 

By this unitary transformation S^^^^^ changes its sign and we can treat the 
N > M > N/2 case. In the case of even N we can show that 

n{-J,-A,h) = U2n{J,A,h)U^\ n 25f = C/2-\ (153) 

l=even 

By this unitary transformation Sf^ Sf , Si change to —Sf, —Sf, Si at I = even. 
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13 Bethe ansatz wave function 



Consider the state where ah spins are up and the total S^' is N/2 

?^|0 >= £;o|0 >, Eo = -JAN/A- Nh. (154) 

Writing a general state \\^ > in terms of a wave function /, as in equation (48), 
the eigenvalue condition can be expressed as 

= -^^(1 -<5„^+i,„,.+i)|/(ni,.,nj- + l,nj+i,.,nM) 

+f{ni,.,nj,nj+i - l,.,n„)| 

+{£;o -E + {JA + 2h)M -JAJ2 }/( ni,n2, .,nM)- 

3 

(155) 

For a wave function of the type eq. (50) to be the eigenstate, set E and A[P) to 
be 

M 

E = Eo + Yyj{^-cos.kj) + h], (156) 

i=i 

= A{P){e'^''= + e-'^^o+i) - 2Z\)e''^o+i) 

+A{P{j,i + l))(e^'=^<^+i) + e-^"-^ - 2A)e^p^, (157) 

A{P) = e{P) ]J(e'('=^'+'=^^) + 1 - 2Ae'^'''). (158) 

l<3 

The periodic boundary condition is as follows, 

AM _ ( 1 ^M-l TT exp[i(fcj- + h)] + 1 - 2^exp(zfc,) 

j = l,2,...,M. (159) 

This is a set of complicated coupled equations for M unknowns. If we have a 
solution of this set, we have one eigenstate and its energy eigenvalue and total 
momentum. If we set rapidity parameters Xj as cot{kj/2) 



exp(zA:,) = f (160) 
\Xj — i) 

the phase factor 



exp[i(A:j + ki)\ + 1 — 2Z\cxp(iA:j) 
exp[i(A;j + A:;)] + 1 — 2Z\exp(ifci) 



Thermodynamical Bethe Ansatz and Condensed Matter 29 

cannot be written as a function of xj — xi except in the case A = 1. This is not 
convenient. If we set 

sin ^(x,- + i) , ^ 

eMikj) = . I) ' (161) 
sm^[Xj — I) 

in place of (160), the phase factor becomes 

cos ^{xj + xi){cosh4> — A) + {Acos^{xj — xi + 2i) — cos ^{xj — Xi)) 
cos ^{xj + xi){cosh.(j) — A) + {Acos^{xi — xj + 2i) — cos ^{xi — Xj)) 

At cosh (fy—A = 0, this phase factor becomes a function of xj —xi and independent 
of Xj + Xi. It is written as 



Then for Zi > 1 (159) becomes 



'sin|(,rj + -j-j- siu§(.rj - .r/ + 2i) 



^ sin %{xj-i)^ f^. sin | {xj -xi-2i) 

(j) = cosh"^ A, ^ > 0. (162) 



For Zi < — 1 we set 



sin ^ (xj +i) , ^ 

eMikj) = - . ■ ■ (163) 



Equation (159) becomes 

'sm^{xj + -f-f sm^{xj - xi +2i) 



Vsin i-(xA — iW 



^ sin ^ {xj — i)' f'_^ sin ^ (xj — xi — 2i) ' 

= cosh-^(-Zi), (/)>0. (164) 

In the case 1 > Z\ > — 1 we set 

, , , sinh ^(x,- +i) , ^ 

exp zfc, = - . J, ^, 165 

smh ^ (Xj —t) 

yielding the Bethe ansatz equations 

/ sinh ^ [xj + „ sinh ^ [xj - xi + 2i) 

Vsinh ^(xj — i) / j-J: sinh ^(xj — x; — 2i) ' 

7 = cos"^(-zi), 7r>7>0. (166) 
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14 String solutions at ^ > i 

By the transformation (161) the wave function and eigenvalue are written as 
follows, 

/(ni,n2, ...,nM) 

V 7<; 2 ^t^^ slnf(a;p,■-^) 

i; = i^o + E(2/.+ /r^'^ ), i^ = E2cot-^^^^fe^- (168) 
cosh 0- cos 0a;/' ^ tanh((?i/2) ^ ^ 

The following string solutions are possible for complex xj 's for the N = oo case 
from the normalizability of the wave function, 

Xj =a + {M + l- 2j)i, (169) 

where a is a real number at —Q < a < Q, Q = 'k/4>- The total momentum and 
energy is given by 

1 tan(0a/2) „ ^ J sinh sinh M(/) ,,„„x 

ii^ = 2cot-^ , E = Eo + — — f-+2Mh. (170) 

tanh(M(/)/2) ' cosh M0- cos 0a ^ ' 



Thus the dispersion is 

" cosh M(j) — cos K 



E = Eo + 2Mh + Jsmhc, 



sinhM0 



(171) 



This excitation energy gives (71) in the limit Z\ ^ 1. In the limit of large A the 
energy is J A + 2Mh. This is M successive down spins in a sea of up spins. The 
lowest energy state at /i > 1 with M down spins is the M string state given by 
(171) with zero total momentum. So the energy of the ground state for Zi > 1 is 

£; = -:^^-(7V-2M)/i+ Jsinh^tanh^. (172) 

15 Thermodynamic equations for the XXZ model for 
A > 1 

Gaudin derived a set of thermodynamic Bethe ansatz equations at Z\ > 1 
(Gaudin (1971)). The wave function for M down spins in the infinite lattice 
is 

M M 

/(m, n2, hm) = {ziZ2..ZMr' E ^(^) n (11 ^PlT'^'~"' > 

P j=2 l=j 

sin 4(x, , , 

= e'"' = i . I)' . )■ (173) 

smf(xj -i) 
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This corresponds to (73). From the normahzabihty condition of the wave func- 
tion we have 



^(7) 7^ 0, A{P 7^ /) = 0, \ziZ2...zm\ = 1 



M 



I n zi\<l, i = l,2,3,...,M-l. 



These conditions are satisfied only if 



a + {M+l-2j)i, j = l,2,...,M, Q>a>-Q. 



(174) 



(175) 



We can show that 

M 



sinf(a + i(M-2j)) 
sin|(a-zM) 



/ cosh 0(M — 2j) — cos 4>a 



< 1, 



cosh (pn — cos (f>a 

1 < J < M - 1, 



for arbitrary M. Thus a string with arbitrary length is possible for the XXZ 
chain at \A\ > 1. In the case |Z\| < 1 the string condition is more complicated 
than the case |zi| > 1. From (162) we have the following equation corresponding 
to (80), 



Here 



sin ^{x + in) 



Enm ('^) 



en(a;) — . , , 
sm ^{x — m) 

e|n-TO|(2;)e^„_„l_i_2(a;)e^^_^l_i_4(a;)...e^_i_^_2(a;)e„+TO(a;) 
for m, 
el{x)el{x)...eln_2{x)e2n{x) for n = m. 



(176) 



(177) 



(178) 



a;^ is the real part of a-th string in the strings of length n. The logarithm of 
(176) is 

Nen{x'^)=2nI2+ Onmix'^-xf), (179) 

{m,0)jt{n,a) 



where 
and 

^nm{x) 



^ / N „ 1 / tan ^ ^ ^ X 
0„(a:) = 2tan-i(-— ^)+27r[- 
tanh -f- 



2Q 



^ n— m (x) + 20\n-m\+2{x) + ... + 20n+m-2{x) + On+m{x) 

for n ^ m, 

202{x) +2e4{x) + ... + 2e2n-2{x) +92n{x) for n = m. 



(180) 
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The function 6n{x) is a quasi periodic function which satisfies 

6n{x + 2jQ) = Onix) + 27ri, j = integer. 

We consider the energy of general eigenstates which is given by the set of quan- 
tum numbers {1^}, 

mi:}) = N{-h -■^)+ ^( '"■'^"^^^^ a„(x^) + 2hn), (181) 

n,a 

where 

1 (psinhrKf) 
^ ^ 277 cosh n4> — cos (j)x 
The partition function of the XXZ model is as follows, 

[^/2] (N-IM^h 
^= ^(l + (l-<5^,2M)exp-^ ^X)^exp[-T-i£;({0)]. (183) 

M=0 {/J} 

Corresponding to (88) we define the following functions, 

K{X) = 9n{x) - iV-1 ^nm{x - X^). (184) 

(m,a) 

Using this function we can determine the position of holes for n-strings. We define 

the distribution functions of particles and holes of n-strings as Pn{x) and Pn{x). 
By the equation (179) we have the conditions for these two kinds of functions 

{x) + Pn{x) +YTnm* Pm{x)- (185) 

m 

Here was defined in (182) and 

{^\n-m\{x) + 2a|„_„|_,_2(a;) + 2a|„_„|+4(a;) + ... 
+23in+m-2{x) + an+m{x) for n 7^ m, (186) 
2a2(a;) + 2a4(a;) + ... + 2a2n-2(a;) + a2„(a;) for n = m. 

Here the meaning of convolution of two periodic functions a and b with period- 
icity 2(5 is redefined 



-Q 

The energy per site is 



rQ 

a*b(a;)= / ai{x - y)h{y)dy, (187) 
J-Q 



J A °° f'^ 

e = -(— + /i) + y] / ^n{x)pn{x)da 

4 rrr J-Q 



, , 27rJsinh(/) , . ^ , ,^^r.-. 
gr,{x) = — -an{x) + 2nh. (188) 
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The entropy per site s is 



OO pQ h / \ 

E/ Pn(-)Kl + ^)+p!:(x)ln(l + ^ 

71=1 "'-Q f>n\-i') Pn 



)dx. (189) 



The condition of minimizing the free energy e — Ts yields equations for rin{x) = 

Pn{x)/pn{x), 

ln,7„(x) = ^ +Y.^nm* ln(l + r?-^(x)). (190) 

771=1 

This set of equations is equivalent to the following one 

, , 27rJsinh(i) , . , ,^ , ., , , 

lnr?i(ar) = ^^^^s{x) + s * ln(l + r}2{x)), (191) 

\nr]n{x) = s * ln(l + rjn-i{x)){l + Vn+i{x)), (192) 
Inrjn 2h 

hm — — = — , (193) 

71— >oo n 1 

where 

s(x) = I E sech(!^^^^). (194) 



4 

n=—OG 



The free energy per site is 

.JA 



+ /i)-TV/ ain{x)Hl + v-\x))dx. (195) 



Corresponding to (104) we have another expression for the free energy, 

f = J , ^ / ai ix)s(x)dx T s{x) ln(l + r]i (x))dx. (196) 

(115-117) and (191-193) have the almost same structure. These equations are 
called Gaudin-Takahashi equation (Gaudin (1971), Takahashi (1971a)). 



16 Theory for |^| < 1 XXZ model 

16.1 String solution of an infinite XXZ model with |.^| < 1 

The shapes of strings for \A\ < 1 are quite different from those at > 1. Taka- 
hashi and Suzuki proposed a condition of the strings and constructed thermody- 
namic integral equations (Takahashi and Suzuki (1972)). Later Hida, Fowler and 
Zotos derived these conditions from the normalizability condition of the string 
wave function for an infinite chain (Hida (1981), Fowler and Zotos (1981)). For 
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|Zi| < 1 there are two kinds of strings, one of which has the center on the real 
axis and the other is centered on the poi axis, 

Xj = a+ {n + 1 - 2j)i, j = l,2,...,n, (197) 
Xj =a+{n+l-2j)i+poi, j = l,2,...,n. (198) 

We designate that the string of the former type has parity v = 1 and that the 
latter has parity v = —1. Applying the normalizability condition of the form 
(197) yields 



sinh ^(a + i{n - 2j)) I cosh 7a - cos7(n - 2j) 

i=j+i 



^ ^ n "^'1 ' sinh^(a — in) ' V cosh 70: — cos 771 



Thus cos7n < cos 7(71 — 2j) for j = 1, 2, 3, n — 1. For (198) we have 



^ -pi- cosh ^ (a + i{n — 2 j)) / cosh 70; + cos 7(71 — 2j) 

cosh^(a — m) y cosh7a + cos7n ' 

and therefore cos 771 > cos 7(77 — 2j) for j = 1,2,3, ...,n — 1. Then from the 
normalizability condition we get 

< v{cos{{n — 2,7)7) ^ cos(n7)) 

= 2z;sin((n — j)7) sin(j7), for j = 1, 2, n — 1. (199) 



This is equivalent to 

(_1)[^^]+[^J for j = l,2,..,n-l, (200) 

for j = l,2,...,n-l, (201) 

TT TT 

where [x] denotes the maximum integer less than or equal to x (Gauss' symbol). 
For rational po = 7r/7, (201) is a strong condition. If po = nxjni^ and rii and 
77-2 are coprime, the string with length greater than n\ cannot satisfy at least 
one of (201). Thus n > ni + 1 strings are forbidden. Moreover for a n = n\ 
string, the momentum is always or tt. So this string has also no meaning for 
the thermodynamics. Next we seek the number n and parity v which satisfies 
(200) within n <n\. Equation (200) is equivalent to 

[^^^] + [^] - [ + [i^](Mod2), , = 1, 2, n - 2, 

TT TT TT TT 

. 1^ (Mod2). 

As [i^^] - [ ^""7'^^ ] is or 1 and [^] - \^^^\ is or -1, we obtain 
[(!i^] + [j2] = [ (--i-l)7 ] + [(l±l)7]^ ,- = l,2,...,n-2. 

TT TT TT TT 
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These are strong restrictions on the parity v and the length of the string n. The 
above conditions are eqiiivalcnt to the following conditions which were given in 
(Takahashi and Suzuki (1972)). The length n of a string should satisfy 

n-l 

2 J2 = ('^ - 1) - ^hM ' (202) 

t'sin{(n- 1)7} > 0. (203) 

This condition was first introduced under the assumption that these strings 
form a complete half-filled state (Takahashi and Suzuki (1972)). Later Hida, 
Fowler and Zotos showed that conditions (202,203) can be rcdcrivcd from the 
nor maliz ability conditions of the wave function for N —>■ 00 and finite M (Hida 
(1981), Fowler and Zotos (1981)). For a given value of A (or 7) we can determine 
the scries of n which satisfies the conditions (202) and (203). If 7/7r is a rational 
number, this series becomes finite and the number of unknown functions is also 
finite. We consider the 7 = w/u, v = integer case. In this case n = l,2,...,j/ 
satisfy (202). For n = 1 both v = 1 and —1 arc possible. For n = 2, 3, — 1, 
only V = +1 states are possible. These excitations have the following energy and 
momentum 

E = -2 J + 2nh, (204) 

t;cosh(a7) — cosn7 

j^=-nn(-""^fi"^ + ^;^-^^"{' + ^"^;). (205) 
V sinh ^(q!7 + i(l — w)7r/2 — zn7) / 

The energy and momentum have the following relation, 

„ ^ . cos 717 — cos K „ , 
E = -Jsmj \-2nh. 

smn7 

The momentum is restricted to the region 

7^7 

\K\ < TT — {nj — n[ — ]) for v = 1, 

717 

> \K\ > TT — (717 — 7r[ — ]) for V = —1. 

TT 

Then for n = v, the energy and momentum are always zero. Only one state 
is obtained from this string solution. So we exclude this n = v state from the 
thermodynamics of this case. So v string states (1, +), (2, +), {f — l, +), (1, — ) 
play important roles. Especially at Z\ = 0, 7 = 7r/2, u = 2 we have only string 
states (1,+),(1,— ). These are single states at momentum \K\ < n/2 and \K\ > 
■k/2. 

Next we consider the 7 = Tr/{vi + 1/1^2) case. (1, +), (2, +),..., (i^i — 1,+), 
(1, -), (1 + /.!,+). (1 + 2i.i, -), ...(1 + {V2 - (-1)"^-'), {i^u {-lY') satisfy 
conditions (202, 203). Thus 1^1 + V2 strings are necessary to describe the ther- 
modynamics of this case. 
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For a general rational number between and 1, we can express it by a 
continued fraction with length I, 

- = ^ + ^ + + ui,iy2,...,ui-i>l, ui> 2. (206) 

We define the following series of numbers y-i,yo,yi, ■■■,yi and 
mo, mi, ...,mi as 

y-i=0, 1/0 = 1, 2/1 = 2^1 and yi = yi-2 + i^iVi-i, 

i 

TOo = 0, mi = ^Vk- (207) 
fe=i 

The general rule to determine the parity v and length n is as follows 

rij = yi-i + {j - m,)y„ Vj = (-1)[K-i)/po] for mi < j < mi+i, 

rim, =yi-i, Vm, = (-1)'. (208) 

The number of strings is mi. We give examples for some rational numbers in 
Tables (1), (2). We put x" as the real part of strings with parity vj and length 
rij . a takes values from 1 to Mj . We find the following relations for these series 
of numbers 

= ^[(1 ^ 2(5„i^ j)nj_i + rij+i], for < j < mi+i - 2, 
Tij = (1 - 2dm,-ij)nj-i + rij+i, for j = m^ - 1, i <l, 
no = 0, rimi + rimi-i = yi- (209) 



Table 1. Length rij, parity Vj and Qj of strings for some rational values of 'y/n 



j 


i 

6 


3 
16 




1 


13 
69 




i 


5+1 




1 


1 


+ 


4 


1 


+ 


13/3 


1 


+ 


56/13 


2 


2 


+ 


3 


2 


+ 


10/3 


2 


+ 


43/13 


3 


3 


+ 


2 


3 


+ 


7/3 


3 


+ 


30/13 


4 


4 


+ 


1 


4 


+ 


4/3 


4 


+ 


17/13 


5 


1 




-1 


1 




-3/3 


1 




-13/13 


6 








6 


+ 


-2/3 


6 


+ 


-9/13 


7 








11 




-1/3 


11 




-5/13 


8 








5 


+ 


1/3 


5 


+ 


4/13 


9 














21 




3/13 


10 














37 


+ 


2/13 


11 














53 




1/13 


12 














16 




-1/13 



Thermodynamical Bethe Ansatz and Condensed Matter 37 



Table 2. Length rij, parity Vj and qj of strings for conjugate values of y/ir in previous 
tabic 



j 


5 


— 




1^ 

16 


— 




69 


— 




i+I 








4-1 


1 


1 




-4/4 


1 




-13/13 


1 




-56/56 


2 


2 


+ 


-3/4 


2 


+ 


-10/13 


2 


+ 


-43/56 


3 


3 




-2/4 


3 




-7/13 


3 




-30/56 


4 


4 


+ 


-1/4 


4 


+ 


-4/13 


4 


+ 


-17/56 


5 


1 


+ 


1/4 


1 


+ 


3/13 


1 


+ 


13/56 


6 








6 


+ 


2/13 


6 


+ 


9/56 


7 








11 


+ 


1/13 


11 


+ 


5/56 


8 








5 




-1/13 


5 




-4/56 


9 














21 


+ 


-3/56 


10 














37 




-2/56 


11 














53 


+ 


-1/56 


12 














16 


+ 


1/56 



16.2 Scattering pheise shift among strings 

Corresponding to (80) and (176) we have Bethe ansatz equations for strings as 
follows 

nil Mk 

= - n n - (210) 

fe=l/3=l 



where 



ej(a;) = 9{x;nj,Vj), 



(211) 



5(a;;2nj,'yjUfc)nr=i 9'^{x;2l,v.jVk) for rij = Uk, 

,1 
) 

5(x;n,+) = "--^- ' -^--^ ^- 



Ej,k{x) = { g{x\ (rij +nk),VjVk)g{x; jn^ - nk\,VjVk) 

X g,2(a;; \nj - nfcl + 21, vjVk) for nj n^, 

sinh ^{x + in) cosh ^ (a; + in) 

g{x;n,-) = -- 



sinh ^{x — in) 
The logarithm of (210) yields 



cosh ^{x — in) ' 



mi Mfc 



where 



fe=l /3=1 



= fix;nj,Vj), Ojkix) = f{x; {n^ - nk\,VjVk) + 

Min{nj — 1 

f{x;nj +nk,VjVk) + 2 ^ f{x;\nj - nk\ + 2i,VjVk), 

i=l 



(212) 
(213) 

(214) 



(215) 
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a) 7 = 7r/5 



b) 7 = 47r/5 



c) 7 = 37r/16 



d) 7 = 137r/16 



16 



16 
13 



e) 7 = 137r/69 



f) 7 = 567r/69 



69 
13 



69 

56 



Fig. 1. Strings on the complex plane for given values of po = 7r/7. If we change 
7 TT — 7, the strings are almost the same. But imaginary parts of -yxj shift by 

TT. 
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and 

f(x- n v)^l ° ^"""^ ""Z^'" = (216) 

jyx, ri, u) 2utan-i{(cot(n7r/2po))" tanh(7rx/2po)} otherwise. ^ ' 

The quantity I,-^ is an integer (half-odd integer) for Mj odd (even), which is 
located in the region 

-. mi 

< ^WOj{oo) - ^Mfe0,-,,(oo)|. (217) 

The function f(x:nj,Vj) is a monotonically increasing function for m2i < i < 
7712 j+i and a monotonically decreasing function for TO2j-i < j < TO2i. 

16.3 Bethe-ansatz equation for XXZ model with |^| < J 

Following Yang and Yang we define particles and holes of strings. We obtain an 
integral equation for distribution functions pj and p'j of particles and holes of 
strings in the thermodynamic limit, 

mi 

ajix) = sign{qj){pj{x) + p'jix)) + ^ T,-,fe * pk{x). (218) 

fe=i 

Here 



aj{x) s {27r)-'-0j{x), Tj,k{x) ^ (27r)-i-0,- ^(a;). (219) 



-lo.ix), T,,(.).(2.)-i 
The symbol a*b denotes the convolution of a{x) and b{x) as follows 



a * b{x) = a{x- y)b{y)dy. (220) 

J —CO 

The functions aj{x) and their Fourier transforms aj{u)) are written as 

/ N 1 7 sin 70," , , , sinhowo; 

aJx) = '-^ , a, w = . , ^ , 221) 

27r cosh 7.T + cos 7(7j smnpo"-^ 

Qj = i-iy{pt - (j - mi)pi+i), for mi<j< nii+i, (222) 
where the p,, i = 0, 1, I are defined by 

PO = 71-/7, Pi = 1, ^^i = bi-lM]> Pi =Pi-2 -Pi-lt-i-l. (223) 

For the series qj we find the following relations 

9i = ^[(1 - 2(5m,j)'7j-i + for < j < m^+i - 2, 

«j = (1 - 2(5to,_i j)gj_i + for j ^rrii-l, i < I, 

qo = Po, rimi + n-mi-i = 0. (224) 
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The Fourier transform of Tj^h (x) is given by 

o- / ^ t-ut . sinh((po - IgjlV)sinh(g'few) 

= 2sign gj)coth(p,+iw) ■ ' ; \ ^ 

sinn(pot^j 

+Sj,mi-i5k,mi - Sj,k, for j <k, rrii < j < rui+i. 
At j = 1 we have 

Ti,k{^^) = sign(gi)2cosha;afc(a;) - ^i,fe. 

At j = mi we have 

The energy and entropy per site are given by 

,JA '•^ 



+ h) + Y^ / 9j{x)Pj{x)d2 



, . 27rJsin7 
gj{x) = o,3\x) + zrij/i, 



and 



Pj ln(l + ^) + p'] ln(l + ^)dx. 

-oo pj Pj 

To minimize the free energy density e — Ts with respect to pj , we have 
S{e -Ts) = Y^ r gAx)5pAx) - tI^p,- ln(l + ^) + ln(l + ^) 

The variation of (218) gives 

= - sign(gj) ^Tj-fe * Spk- 

Thus we obtain 

/oo 
-oo 



h 



{M^ _ in(^) + ^sign(g,)T,,, * ln(l + ^)}. 
^ Pj k 



Thermodynamical Bethe Ansatz and Condensed Matter 41 



At the thermodynamic equilibrium one obtains the following non-linear equa- 



tions determining r]j{x) = pj{x)/pj{x), 



lmjj{x) = gj{x)/T + ^sigu{qk)Tk J *ln{l + T]^\x)), j = l,...,mj. (230) 



fe=i 

The free energy is given as follows, 



f = e-Ts = -{^+h) 



+ E / PiWtffi W -Tlnvjix)] - T[pj + p"^] ln(l + 777!). 

If we substitute (230) into the first bracket and (218) into the second, the Tj^^ 
terms are cancelled and we get 



/ = -(— + /i)-T^sign(g,) / aj{x)\n{l + 7ij\x))dx. 
If one uses the j = 1 case of equation (230) and (226), one obtains 



(231) 



. , , 27rJsin7 p , , , 

f = — 7 sign(gi) / ai{x)si{x)dx 

* T J —00 

/oo 
si{x)\ii{l + r]i{x))dx. (232) 
-00 

Prom equation (221-223) we get the following relations, 

aj - Si* ((1 - 2(5m,_i j)aj_i + %+i) = 
for TOj-i < j < rui - 2, 

arrn-i - (1 - 2(5mj_i,TO4_i)si * ami-2 - di* arm-i - Si+i * drm = 
for i < I, 

ami-i{x) = -Umiix) = Si * ami-2, (233) 

where 

ao(a;) = S{x), 

f°° duj e^"^^ 1 , -Kx 

sAx) = / — — - — sech- — , 

J_oo 47r cosh(pia;) Api 2pi 

c^.(.)^ r^ ^^7;^'f T^^^^t (234) 
7_oo 47r cosh(piw) cosh(pi+iw) 
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Using (225) one can show the following relations, 
Tj,h - Si* ((1 - 25to._j + Tj+i^k) 

for rrii-i < j < — 2, 

— Si+1 * Tmi,k = { — ^y^^ {Smi-2,kSi + Smi-l,kdi — Smi,kSi+l), 

for i = 1, 2, ., Z — 1, 

Tmi-l,k = -Tmi,k = Si* Tmi-2,k + sigll{qk)Smi-2,kSl, (235) 

with To,k = 0. Using (233) and (235) one can rewrite (218) as follows 
Pj +p'} = si* (p^_i + p^+i) for rrii-i < j < - 2, 
Prrn-l + pti-i = Si * pti-2 + di* pti-i - Si+i * pti, 

Pmi-l + p'^i-1 = Pmi + Pm, =Sl* (236) 

with Pq = 6{x). Equations (230) are rewritten as 

ln(l + .o) = -?^^(x), 

7-1- 

lnr]j = (1 - 26m,.i,j)si * ln(l + ijj^i) + Si* ln(l + r]j+i), 

for rui-i <j<mi-2,jy^mi-2 

lnr?TO,_i = (1 - 26mi_-^,mi-i)si * ln(l + r?m,_2) 

+di * ln(l + r?m,_i) + Sj+i * ln(l + r;„J, for i<l 

lnrimi-2 = (1 - 25m,_i,m,-2)si * ln(l + r]m,-3) 

+s,*ln((l + r7„,_i)(l + r,-;)), 

lnr;m,_i - yih/T = yih/T - In rjmi 

= s;*ln(l + r?„,_2). (237) 

Then if we write ln«;(a;) = \nr]mi-i — Vih/T we have integral equations with 
mi — 1 unknown functions 

ln(l + ,o) = -^^^(.), 

ln?7j = (1 - 25m,.i,j)si * ln(l + r]j-i) + Si * ln(l + 77^+1), 

for rrii-i < j <mi -2,j ^ mi - 2 

lnr?m,_i = (1 - 25mi_umi-i)si * ln(l + r?mi-2) 

* In(l + ?7m.-i) + Sj+i * ln(l + 77mJ, for i < Z 

ln/7m,-2 = (1 - 25mi_i,mi-2)si * ln(l + 7?TO,-3) 

+s; * ln(l + 2 cosh{yih/T)K + k^), 

lnK(a;) = sj * ln(l + r7m,_2). (238) 
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17 Some special limits 

17.1 T — > cx> or J ^ limit 

In equations (238), ln(l + r]o) becomes zero and rij{x) are all independent of x. 
This yields the following difference equation 

v^ = {i + vj-i)'-''-^-^Hi + Vj+i) 

for TUi-i < j <mi-2,j ^ mi -2 

Vl,-i = (1 + r?m,-2)'"'^'"*-^''"*-' (1 + Vm,-i) ln(l + 
for i < I 

<-2 = (1 + '7m,-3)'"'''"'-— ^(1 + 2 COSh{yih/T)K + K^), 

= (l + ??m,-2). (239) 

The solution of this set of equations is 

_ sinh(nj + j/^-i)fe/T 2 ^ 
^ sinh(yi_i/i/r) ' 
for rrii-i < j < rrii, j <mi- 2, 
^ sinh(n,n,-2 + y;)/i/r 
sinh(y,/i/T) 



(240) 



For j = 1 wc have r?i = {2coshh/Tf - 1. Substituting this into (232) we find 

the free energy 

f/T = - l-a{2coshh/T). (241) 
At /i = this gives that entropy per site is In 2, as it should be. 



17.2 J > 0, T ^ limit 

We define ej{x) = T\nrjj{x) and e^(a;) = Tln(l + rjj{x)). One can show that 
Ejj i > 2 is always positive. The equation (230) gives 

27rJsinT 1"^ 
ei(a;) = '-a{x,l) + 2h- a{x - y,2)e:[{y)dy. (242) 

T J — oo 

If e{x) < at |a;| < iJ and e{x) > at > B, then one obtains a hnear integral 
equation for pi{x), 

pi{x)+ j a{x-y,2)pi{y)dy = a{x,l), (243) 

J-B 

where 

1 7sinn7 

a(x,n) = - . (244) 

27r cosh jx — cos 



5. Thermodynamics of the XYZ model 



18 Bethe-ansatz equations for the XYZ model. 

Here we consider the symmetry of the following Hamiltonian 

N 

H = - ^ JxSiSi^i + JySfSf_^-^ + JxSfSij^-^. (245) 

1=1 

We assume N is even. By the transformation 

l=even 

Ti^Jx, Jy, Jz) 'H{—Jx, ^Jy, Jz)- In the same way 

W(Jx, Jy, Jz) —>■ Ti-iJx, -Jy, -Jz) and H{Jx, Jy, Jz) Ti-iJx, -Jy^ -Jz)- Namely 
the energy spectrum of this Hamiltonian is unchanged for reversing signs of two 
Jq,'s. It is evident that the spectrum is unchanged for exchanging J^'s. Thus it 
is sufficient to treat only the case 1 > Jy/Jz > \ Jx\/Jz > 0. Baxter solved the 
eight-vertex model and also the XYZ model (Baxter (1972a), Baxter (1972b)). 
The Bethe ansatz equation for this model is 



N/2 

, HiiiaXi+i)) \N ^ 2niu'/po TT HiiiaXi-Xj + 2i)) 

^ Hi {iCixi -i))J Hi {iC{xi - Xj - 2i)) ' 

N/2 

Y,xi=Qv' + ipov, Q = K{l')/C, po = K{l)/C. (246) 
Here the modulus I and the parameter ( are determined by 



' = yi|Tri' cn{2c,i) = -Jx/Jz. 

There are N/2 rapidities. Hi{x) is the Jacobian elliptic function defined by 
Hi{x) ^ 2 f;(-l)"+ig"("-i)+V4 3in(2n - 1)^, exp(-^!|^). 
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This function has the foUowing properties, 

Hi{x) = -Hi{x + 2K{1)) = -ge*"^/^(')i?,(x + 2iK{l')). (247) 

In the Bethe ansatz equation (166), the function sinh^a; is merely replaced by 
the elliptic theta function. The energy is given by 

N/2 

E = -^[l- ^(a(0, 1) + a(g, 1))] - '-^Y.<-u D, (248) 
where 

^ ' 2mdx \Hi{iC{x-il))J ^ ' 

19 Strings and the thermodynamic Bethe ansatz equation 
for the XYZ model 

In the limit I equation (246) becomes 

7V/2 

. sm(/c,(.]7 +^)) \^ _ ^-2^iv'lv„ TT «m(^(-' f - j + ^0) 
^sin(iC(a;i - i)) J sin(zC(a;; - xj - 2i)) ' ^ ' 

This equation is equivalent to (166), if we assume that the x;'s are finite and 
u' = 0. ( becomes 7/2 and Ki becomes n/2. So it is natural to assume the same 
types of strings can be determined using po = Ki /(, 

Xj = a + (n + 1 — 2j)i, 

Xj = a + {n + 1 - 2j)i + poi, Q>a> -Q. (251) 

The solution becomes doubly periodic. So we consider the distribution of solu- 
tions at —Q < ^x < Q and —po < "^x < po- It is expected that the same kind 
of strings appear in the case of XXZ model at 7r/7 = po- We can determine rij's 
and qj's via (209) and (222), 

m, Ah 

{e,(4)}^ = -exp(-27rir.7po) J] 11 ^^'fc^^" - 4)' 

k=l 0=1 

^'=nEE".-^'- (252) 

^ j=l a=l 
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where 



ej{x) = g{x;nj,Vj), (253) 
Hi{i({x + in)) 



g{x;n,+) - . 

Hi{iQ[x - m)) 

Hi{Ki + iC{x + in)) ^254) 




HiiKi+iCix-in))' 

g{x;2n.j,VjVk)Y[ili^ g'^{x;2l,VjVk) for rij = nk, 

g{x] {rij + Uk), VjVk)g{x; \nj - nk\,VjVk) , . 

X n;^f g\x; \nj -nk\ + 21, v,Vk) ^ ' 

for rij 7^ nfe, 

Taking the fogarithm of (252) we have 

N9j{xi) = 2nli-2TTiy'/po + ^^Ojk{xi-x'l), a = 1,2,..., Mj. (256) 

fe=l/3=l 

Here 

9j{x) = f{x; nj,Vj), Ojk{x) = f{x; \nj - nk\, vjVk) + 

Min(nj .n^) — 1 

f{x;nj+nk,VjVk)+2 ^ f{x;\nj - nk\ + 2i,VjVk), (257) 
and f (a;, n, v) is defined by 

CxD 

f (a;, n, v) = f{x, n, v) + f{x — 21Q, n, v) + f{x + 21Q, n, v). 

1=1 

f{x,n,v) was defined in (216). An eigenstate should be identified by the set of 
quantum numbers I". Prom (248) the energy must be 

T rv/- T ™i ^-fj 

^ = -ATJ^J? — ^ ^a(a;i, 1) = -A^ J,i? - 2-1 

^ i=l ^ j=l a=l 

i«=J[l-^(a(0,l)+a(Q,l))], 

a. (-) - + 2 E ^^^^^^ cos(.j./Q)l , (258) 

2QLpo ^ sinh(po7r//<3) ^-^'^^J' ^ ^ 

where qj was defined in (222) The number of zeros must be N/2, so 

mi 

N/2 = J2njMj. (259) 
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Thus the energy per site is given by 



e = -J^R - 
The entropy per site is 



^.,{x)p,{x)dx. (260) 



s = Y, Pi 1^(1 + -) + 111(1 + (261) 
Prom (257) we have the relation between pj{x) and Pj{x), 

mi 

aj{x) = sign{qj){pj{x) + p){x)) + ^ T,- fe * pk{x). (262) 

fe=i 

Moreover, from (259), 

- = m = ^ nj y p, {x)dx. (263) 



Next we need a Lagrange multiplier to guarantee the condition (263). One should 
minimize e — Ts + 2hm under conditions (262), and after that the multiplier h 
should be chosen so that (263) is satisfied. Just in the same way as before we 
get the integral equations for r/j(x) = pj{x)/ pj{x), 

mi 

lurjjix) = gj{x)/T + Y^sign{qk)Tk,j *H1 + r]^:\x)), j = I,..., mi. (264) 

k=l 

Here *, gj(a;) and Tj^k{x) are 

f * g{x) = f{x- y)g{y)dy, 
J-Q 

gj{x) = -^^^^^^^a.j{x) + 2njh, 
T,-,fc(x) = — Yl e''^""/'^f,. fe(^) = J2 TjMx-21Q). 

1 — — 00 1 — — 00 

The quantity g = e — Ts + 2hm is given as follows 
g{J,,T, h) = -J,R 

+ E / pA^) Sj{x) - Tlnr]j{x) -T pj+ p) ln(l + ^T^)dx 
^^^J-Q L J L J 

nil ,.Q 

= -J,R - T ^ sign(g,) / aLj{x) ln(l + r]J^{x))dx. (265) 

A — ^ J —O. 



48 Minoru TAKAHASHI 



Corresponding to (232) this is 

g = -J2/?+ /i - sign(gi) / ai(a;)si(a;)da; 

^ J-Q 

rQ 

-T / si{x)ln{l + r]i{x))dx. 

J-Q 



(266) 



Then m should be determindcd by 

The equation (264) is also equivalent to the following block tridiagonal equations, 

, N 7rJ2sn2C , 
ln(l+r?o) = ^S{x), 

Inrjj = (1 - 2Smi-i,j)si * ln(l + r]j_i) + Sj * ln(l + rjj+i) 
for mi_i <j<mi-2,j^mi- 2, 

ln77™^_i = (1 - 2(5„,_j_m,_i)si * ln(l + r7r„,-2) 
+di * ln(l + rjmi-i) + Si+i * ln(l + ?7„J for i < I, 
\nrimi-2 = (1 - 2(5m,_i,m,-2)Si * ln(l + rimi-2.) 

+si * ln(l + 2cosh(?/;/i/T)K + k^), 

InK = s; *ln(l + r?„,_2). (268) 

In this equation parameter h appears only in cosh{yih/T) term. So rii{x,h) is 
even function of ft- and '^^gj'^^ \h=o — 0. Using (267) we find m = 1/2 and the 
condition (263) is satisfied at ft = 0. Thus equation (268) becomes 

ln(l + ?7o) = |^(5(a;), 

ln?7j = (1 - 26m,.i,j)si * ln(l + 7]j_i) + * ln(l + 77^+1) 
for mj_i <j<mi-2,jj^mi- 2, 

lnT?mi-l = (1 - '^Smi_i,mi-l)Si * ln(l + Vmi-2) 

+dj *ln(l + 7?„,_i) + Si+1 * ln(l + ?7„J for i < I, 

ln?7m,-2 = (1 - 2(5m,_i,m,_2)s; * ln(l + JJmi-s) 

+2si *ln(l + k), 

luK = s; *ln(l + ?7„,_2). (269) 
Corresponding to (232) the free energy is 

/ = - JzR - sign(g'i) / ai(a;)si(a;)c<a; 

^ J-Q 
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rQ 

-T / si(a;)ln(l +r/i(a;))da;. 
J-Q 

(270) 

We can calculate free energy of XYZ model in zero external field. 



6. Numerical calculation and recent 
developments 



In this lecture we restricted ourselves to the unnested Bethe ansatz. Fermions 
with 5-function interactions and the Hubbard model belong to the nested Bcthe 
ansatz which has several kinds of rapidities. For details see (Takahashi (1971b), 
Lai (1971), Takahashi (1972)). But the method of derivation is essentially the 
same as the unnested cases. 

The equations introduced in this lecture note have been solved numerically. 
One can calculate the specific heat, magnetic susceptibility, magnetization curve 
for the XXZ model, and the specific heat for the XYZ model at po = 2,3 (Taka- 
hashi (1974a)). If we increase po the problem approaches to the A = —1 case. 
For the analysis of thermodynamic quantities not at low temperature, wc can 
use the high-temperature expansion method or exact diagonalization. For the 
investigation of low-temperature thermodynamics of solvable models the Bethe 
ansatz method is the only way. For the spin 1/2 ferromagnetic XXX chain, the 
susceptibility diverges as T~'^ and the specific heat behaves as T~". The esti- 
mations of exponent 7 had been done by many authors. Baker et. al. estimated 
7 = 1.66 using high temperature series expansions (Baker et al (1964)). Lyklema 
obtained 7 = 1.75 using quantum Monte Carlo calculations (Lyklema (1983)). 
By the numerical calculation of thermodynamic Bethe ansatz equations it was 
established that 7 = 2 and a = 1/2 (Takahashi and Yamada (1985), Yamada and 
Takahashi (1986), Schlottmann (1985)). This investigation continued to the spin- 
wave theory for low-dimensional magnets (Takahashi (1986), Takahashi (1987), 
Takahashi (1989)). 

In many cases of one-dimensional quantum systems, one can define the quan- 
tum transfer matrix. The largest eigenvalue of this matrix gives the free energy. 
The ratio of the largest and the second largest eigenvalues gives the correlation 
length. This method was developed mainly by Japanese theorists within the last 
ten years. The XYZ model and the Hubbard model were investigated by this 
method, for other solvable models this method is expected to be applicable. 
Logarithmic anomalies in the low-temperature susceptibility of the XXX anti- 
ferromagnet were found by the numerical calculation of Bcthe ansatz equations 
for the largest eigenvalue of the quantum transfer matrix (Eggert, Affleck and 
Takahashi (1994)). 
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